Origin of 'end of aging' and sub-aging scaling behavior in glassy dynamics 
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Linear response functions of aging systems are routinely interpreted using the scaling variable 
iobsAw, where t w is the time at which the field conjugated to the response is turned on or off, and 
where t b s is the 'observation' time elapsed from the field change. The response curve obtained for 
different values of t w are usually collapsed using values of /i slightly below one, a scaling behavior 
generally known as sub-aging. Recent spin glass Thermoremanent Magnetization experiments have 
shown that the value of jj, is strongly affected by the form of the initial cooling protocol (Rodriguez 
et al., Phys. Rev. Lett. 91, 037203, 2003), and even more importantly, (Kenning et al., Phys. Rev. 
Lett. 97, 057201, 2006) that the i w dependence of the response curves vanishes altogether in the 
limit i bs 3> tw The latter result shows that the widely used i bs/tw scaling of linear response data 
cannot be generally valid, and casts some doubt on the theoretical significance of the exponent /x. 
In this work, a common mechanism is proposed for the origin of both sub-aging and end of aging 
behavior in glassy dynamics. The mechanism combines real and configuration space properties of 
the state produced by the initial thermal quench which initiates the aging process. 

PACS numbers: 65.60.+a, 05.40.-a,61.43.Fs 



I. INTRODUCTION 

In glassy systems, a thermal quench initiates a so called 
aging process, whereby physical observables, e.g. the en- 
ergy, slowly change as a function of 'age', the term con- 
ventionally denoting the time t elapsed from the quench. 
In non-stationary processes, conjugated linear response 
and correlations functions generally depend on two time 
arguments, e.g. the system age t, and its value f w at the 
moment where the field is switched on, or off. However, 
both functions appear to actually depend on a single scal- 
ing variable, namely t ohs /t^, where t ohs = t-t w [H0,||. 
The 'observation' time t t> s , (widely denoted by the sym- 
bol t in the literature) is often used in lieu of the sys- 
tem age as an independent time variable. The expo- 
nent \x is generally near one, and the terms super-aging, 
sub- aging and pure or full aging are used to discrimi- 
nate between cases with /j, > 1, jj, < 1 and jj, = 1, re- 
spectively. Sub-aging (henceforth SA) has mainly been 
observed in spin-glass Thermoremanent Magnetization 
(TRM) data @, |3,M@|, usually after subtracting a 'sta- 
tionary' term which describes the response immediately 
after the field is cut. 

As a quench is unavoidably carried out at a finite cool- 
ing rate, the point t = on the time axis eludes a sharp 
definition in an experimental setting. Possibly as a con- 
sequence thereof, the scaling form of response functions 
is, in spite of a protracted debate 0, [H, 0, H, H, [ToL fll|. 
only partially understood. E.g. the physical significance 
of [i (and of the additional time scale it brings along) 
is unclear, even more so in the light of the recent dis- 
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covery that the TRM loses its t w dependence in the limit 
tobs/t w — ► oo [lH. This observation falsifies the long- held 
hypothesis that SA could be a scaling form generally ap- 
plicable to linear response functions. Coined in Ref. 
the term 'end of aging' (EOA) is also adopted here. It 
describes that, for large values of t/t w , the t w dependence 
of the TRM disappears and is replaced a simple logarith- 
mic time decay. The term should not be construed as 
contradicting the fundamental unity of aging dynamics 
which subtends our description. 

In this work, an expression for the TRM decay, Eq. (jHJl , 
is derived which, depending on the value of the model pa- 
rameter x, interpolates between pure aging (PA), SA and 
EOA behavior. Our derivation combines well known con- 
figuration and real space properties of spin glasses and 
supplements them with a model assumption regarding 
the spatial heterogeneity of the system configuration at 
the beginning of the isothermal aging process. Our an- 
alytical expression for the TRM decay is derived by av- 
eraging the magnetic response of independent domains 
over a suitable distribution characterizing their initial 
state. Even though the ratio iobs/^w nowhere appears 
in the treatment, the TRM decay curves produced by 
our Eq. © can, in accordance with standard practice, 
be empirically collapsed in the SA regime using i bs/^w 
as a scaling variable. 



II. LINEAR RESPONSE AND DOMAIN 
HETEROGENEITY 

Thermalized domains whose linear size grows with the 
thermal correlation length [l2| are ubiquitous in aging 
systems with short range interactions. Real space scal- 
ing descriptions [l3[ rely on their properties to account 
for several dynamical properties, e.g. so-called 'chaos' 
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effects. In a physics context, hierarchies were brought 
into focus by the Parisi equilibrium solution of a mean 
field model |14[ , a solution which has ins pire d hierarchi- 
cal models of spin-glass dynamics, e.g. [la ]. However, 
the hierarchical nature of relaxation dynamics in com- 
plex system is a generic feature [la ] which complements 
real space descriptions and which has e. g. been con- 
firmed numerically by solving master equations in fully 
enumerated sets of states [l?], EH • 

Hierarchical models, see e.g. Refs. flijl l20l l2lj . as- 
sume that nested ergodic components [22J exist in con- 
figuration space. In simple cases, the solution of the 
pertinent master equations shows PA behavior [2(|. For 
thermally activated dynamics, this result can be heuris- 
tically explained as follows: Any ergodic component of 
the hierarchy is at any time t indexed by a real valued 
'dynamical barrier' bit). On the (Arrhenius) time scale 
tb(t) = Cexp(6(t)/T) the component is near a state of 
internal thermal equilibrium. The constant C is the fluc- 
tuation time, i. e. the smallest relevant time scale of the 
dynamics. Consider now a system starting in a compo- 
nent with vanishing initial barrier bit = 0) » 0. After 
aging isothermally for a time t — t w , the component 
is characterized by a barrier b(t w ) ex T\nit w / C). The 
dynamics on scales t ohB < t w has character of quasi- 
equilibrium fluctuations within the same component, 
while off-equilibrium processes involving larger compo- 
nents take over for t a b s > t w . The age t w hence sep- 
arates the two dynamical regimes observed for aging 
systems, and constitutes the dominant time scale for 
t > t w . From this observation, pure aging heuristically 
follows. A mathematically more precise route leading 
to the same conclusion relies on the fact that, within a 
hierarchy indexed by energy barriers, only thermal en- 
ergy fluctuations of record magnitude are able to trig- 
ger irreversible changes of ergodic component, or quakes. 
Since the statistics of record-sized fluctuations in a sta- 
tionary series is known analytically [23|, |24| . assuming 
that all physical changes are statistically subordinated to 
the quakes leads to record dynamics (24 , l25l . |26| and to 
analytical formulas, such as Eq. [3j for one and two-point 
averages of aging processes. In spite of their simplicity, 
the above ideas rationalize a large amount of experimen- 
tal evidence, [13, H3, HH . Yet, they neither account for 
SA nor for EOA scaling behaviors. In order to do so, 
the spatial and temporal heterogeneity of sp in g lasses, 
experimentally demonstrated by Chamberlin [29j . must 
be properly taken into account. 

In the present model, independently relaxing domains 
of a glassy system are all endowed with the same type of 
hierarchically structured configuration space. Nonethe- 
less, their respective contributions to the overall linear 
response are different (albeit related) functions of time. 
The assumption is that domains find themselves in states 
characterized by different dynamical barriers at the end 
of the initial quench, or equivalently, at the beginning of 
the isothermal aging process. 

The physical mechanism behind the difference is likely 



related to the way in which the cooling process proceeds 
near the glass transition temperature, see e. g. ref. @,[3(J. 
Here, the spatial heterogeneity is heuristically described 
by a distribution of initial barriers P(b) . We checked that 
a flat distribution supported in the interval (0, b* M ) and 
an exponential distribution with average b* M lead to sim- 
ilar behaviors. The existence of a finite first moment of 
Pib) seems to be the crucial feature. Since the exponen- 
tial form leads to simpler closed form expressions, this 
form is chosen for mathematical convenience. 

In summary the initial state of the aging process fea- 
ture domains described by the dynamical barriers distri- 
bution 

P(6) = i- eX p(-6/6U (1) 

°M 

Let t* M be the Arrhenius time associated to b* M . As we 
shall see, the quantity 



' b* M Htyc) 

controls all deviations from pure aging behavior. Note 
that x > 0, that x — > for t* M — > oo, and that x — > oo 
for t* M I C. 

III. PURE AGING APPROXIMATION OF TRM 
DATA 

The additive contribution to the magnetic response of 
a single domain is in the present theory given by Eq. ([5]). 
In this espression the function Mo(i/t w ) describes the 
TRM of a domain initially in a state with vanishing dy- 
namical barrier. 

The functional form chosen for AIq reflects that record 
dynamics is a homogeneous stochastic process in the sin- 
gle 'time' variable log(t/t w ). By standard arguments, 
all moments of the process, including the average re- 
sponse, admit eigenvalue expansions where log(t/t w ) re- 
places time. The generic term in such expansions is pro- 
portional to it/t w ) Xk , where is the /c'th relaxation 
eigenvalue. In practice, the expansion can be cut after 
few terms and, as shown graphically in the Appendix, 
two terms (one term less than in Ref. [Io| ) already pro- 
vide an acceptable parameterization of the TRM decay 
in the PA approximation. 

Summarizing, the pure aging scaling ansatz for the 
TRM can be written as 

2 

M„(*) = Mi +V(z-1)Y,T ( z " k ~ !) > ( 3 ) 
fe=i k 

where z = t/t Wl and where r\ is the Heaviside step func- 
tion, and Mj is the 'initial' value of the TRM, which for 
simplicity is treated as a parameter. According to the 
formula, the TRM remains constant and equal to Mi un- 
til the magnetic field is cut. The notation a^/Afc for the 
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prefactor is chosen to simplify the form of the rate of 
magnetization change 



braic manipulations, one arrives at 



rTRM,o(Mw) = T «fc ( — 
1 k=i \ tw 



(4) 



The (negative) pre-factors and exponents are for com- 
pleteness tabulated in the Appendix. 



IV. ORIGIN OF SUB-AGING AND END OF 
AGING 



Let T denote the isothermal aging temperature. As 
mentioned, C denotes the smallest relevant relaxation 
time, i.e. the time associated to the smallest energy bar- 
rier in the energy landscape of a single domain. Consider 
now a domain characterized by the initial barrier b* , or 
equivalently, by the Arrhenius time t* = C exp(&* /T). If 
t w < t* , the behavior at time t — t w remains controlled 
by the initial barrier b* and the domain's contribution to 
the TRM correspondingly depends on t/t*. Conversely, if 
t w > t* the initial barrier has been surmounted at f w , and 
the scaling variable is hence t/t w . In real space, the size of 
the domain grows as a function of the dynamical barrier 
b(t). E. g. a power-law growth in the time domain fl2| 
corresponds to an exponential growth in b(t). In our case, 
b(t) = max(Tln(£/C), b*), and domain growth is delayed 
up to the Arrhenius time t*, compared to the case b* = 0. 

Returning to the form of the response, the contribution 
of a domain with initial dynamical barrier b* is 

m(t,t w ,b*) = r)(t w -t*)M (t/t w ) + r)(t* -t w )M (t/t*), 

(5) 

where 77 is again the Heaviside function, where Mq is 
given in Eq. ([3]) and where t > t w . The formula em- 
bodies the key feature of hierarchical relaxation without 
reference to any specific model. Secondly, it introduces 
spatial heterogeneity, as b*, or equivalently, the Arrhe- 
nius time t*, is allowed to differ across the domains. At 
t w = t* , the Heaviside function rj switches between the 
two scaling forms available for the magnetic response of 
a single domain. If the barriers of the different domains 
in the system are all initially near zero, only the first 
term contributes, and the pure aging behavior given by 
Mo(t/t w ) goes through at the macroscopic level. In the 
general case, the TRM decay takes the form 



M(t,t w ) = M (t/t w ) / P(b)db+ (6) 

Jo 

Tln(t/C) . poo 

M (-e- b/T )P{b)db + M l / P(b)db, 

Tln(t w /C) k JT\n(t/C) 



where P(b) is the probability density for a domain with 
barrier in the initial state. 

Inserting Eq. (TTJ) into Eq. ©, and using simple alge- 



M(t,t w ) = Mi 



C 



I I^J )M (t/t w ) 



x ft/tvr 



M (z)z x - 1 dz. 



(7) 



Using the parameterization of Mo given in Eq. ([3]), the 
last expression becomes 



M{t,t w ) = Mo{t/t v 

-x 2 



c 



Xi + x 



(8) 



Pure aging. is achieved in the limit x = 00, i.e. when 
all initial barriers are equal to zero with probability one. 
Sub-aging is present for intermediate values of x. Let 
now A2 be the largest of the two decay exponents char- 
acterizing the pure aging regime. If and only if x < — X2, 
the asymptotically dominant contribution to the TRM 

for large t is proportional to (^) x . To bear this out, 
we first re-write our last equation as 



M(t,t w ) =M (t/tw) 

x 2 



t_\ \ - <H 

Cj ^ Xi+x 



(9) 



The pure aging term Mo(t/t w ) decays the fastest and can 
be neglected. In what remains, the t w dependent term 
having the slowest decay is (t/t w ) X ' I+x . In order for the 
EOA behavior to set in, this term must be much smaller 
than one. E.g. a relative deviation of the TRM curve 
from EOA equal to 1/10 is reached at time 



iEOA = (10)™*, 



(10) 



Since 



\x+\i\ 



is very large when x + X2 ~ 0, the model 



predicts that EOA may occur on a time scale which 
diverges very rapidly with t w . This is qualitativel y in 
accord with the experimental observations of Ref. [111 ]. 
Secondly, when the exponent x is numerically small, the 
expansion (t/C)- x = 1 ~ x\n(t/C) + 0((x\n(t/C)) 2 ) is 
applicable, and the TRM decays in a nearly logarithmic 
fashion for a wide range of t, also in accord with the 
experimental findings. 



V. ON t b s /C SCALING 

Equation {5} features a clear sub-aging behavior with- 
out any reference to the widely used scaling variable t/t^. 
Nevertheless, model generated TRM data can be em- 
pirically scaled in the traditional manner for interme- 
diate values of t and t w . This is checked numerically 
(i) by evaluating Eq. with C = 1, tjj = 10 and 
with all other parameters given in Table. HI and (ii) by 



4 



T=0.60T 



T =0.60 T 



0.75 



i 0.7 



0.65 



0.64 
10"' 



0.75 



1 0.7 



0.65 



0.6 L 



0.94 



0.93 



0.92 



0.91 




0.3 0.4 0.5 0.6 0.7 0.8 0.9 

T/T 



FIG. 1: (Color on line) The Thermoremanent Magnetization calculated according to Eq. (JS| is plotted for t w = 50, 100, 1000 
and 10000 (black, blue, green and red) and for T — 0.60T g versus iobs/iw, (left figure), and versus tobs/iw, (middle). The left 
figure shows that sub-aging is present, and the middle figure that the standard scaling procedure collapses the data. Scaling 
plots for other temperatures are also performed (not shown). In each case a /j, value is estimated as the one providing the best 
collapse. In the rightmost figure, the empirical values of fi are plotted versus T/T g , showing a clear maximum at T/T K — 0.83. 
The line is only a guide to the eye. 



scaling the curves obtained as usually done for experi- 
mental data. The left panel in Figure (fT]) only confirms 
that iobsAw scaling is not satisfactory for the parameter 
values utilized. The four curves plotted all pertain to 
T = 0.6T g . The values of t w used are t w = 50, 100, 1000 
and 10000. The exact same data are shown in the second 
panel of the same figure. They are now plotted versus 
tohs/tl^, with the value of /i chosen as the one judged to 
produce the best data collapse. Additional curves (not 
shown) were similarly obtained by evaluating Eq. ([H]) for 
T/Tg = 0.4,0.6,0.9, and T = 0.93. The empirical \x val- 
ues thus obtained gauge how close the TRM is to pure 
aging. These values are plotted versus T/T g in the right 
panel of the same figure (the line is only a guide to the 
eye). Interestingly, fj, versus T peaks at T/T g = 0.83, 
the very temperature where \x is experimentally closest 
to one [5j. 

In summary, Eq. ([5]) fully contains the standard sub- 
aging behavior widely seen in spin glasses. Furthermore, 
it demonstrates that the applicability of tobsA™ scaling 
does not per se endow the exponent \x with physical sig- 
nificance, which is plainly absent in our case. 



plicable to glassy systems with short range interactions. 
These might include quantum spin glasses, whose critical 
behavior has recently been investigated [3l|, [12] > and ir- 
respective of whether a true equilibrium phase transition 
exists dH or not [3l| . 

The distribution of initial dynamical barriers plays a 
pivotal role in the theory. Arguably, its form depends on 
the cooling protocol, e.g. fast cooling could give a distri- 
bution more sharply peaked at zero, and lead to a relax- 
ation scaling form closer to pure aging. The width of the 
initial barrier distribution is expressed by the exponent 
x, which is experimentally accessible as the logarithmic 
slope of the TRM decay for very large values of t/t w , 
i.e. in the dynamical regime where the t w dependence 
of the data is absent. It should therefore be possible to 
empirically study, via x, how the initial cooling protocol 
affects the distribution of initial barriers and, indirectly, 
the subsequent relaxation dynamics. 



VI. SUMMARY AND OUTLOOK 



In this work, the known scaling properties of off- 
equilibrium linear response functions in spin glasses have 
been accounted for by combining two aspects of complex 
dynamics: the hierarchical relaxation of independently 
thermalizing domains, coupled with the spatial hetero- 
geneity of the initial domain configurations, as defined 
by their initial dynamical barriers. These barriers would 
uniformly vanish for a system conforming to pure aging 
scaling behavior. Our analysis relies on generic properties 
of complex dynamics, and should therefore be widely ap- 
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0.40 


-0.0878 


-2.7575 


-0.0289 


-0.1869 


0.60 


-0.0981 


-3.0945 


-0.0462 


-0.2609 


0.83 


-0.1365 


-2.9016 


-0.0527 


-0.3102 


0.90 


-0.1309 


-3.2785 


-0.0418 


-0.3098 


0.95 


-0.1117 


-3.4737 


-0.0296 


-0.3495 



TABLE I: The first column contains the ratios of the isother- 
mal aging temperature T at which the measurements are 
taken to the critical temperature T g . The other columns con- 
tain the pre- factors and exponents of the two power-law terms 
appearing in Eq. @- All values are obtained by fits (not 
shown) of quality similar to Fig, ((2j). 



VII. APPENDIX 

We describe in this Appendix how the parameter val- 
ues entering Eql3] are estimated by fitting to experimen- 
tal TRM data. The data are obtained according to a 
standard procedure: a Cuo.94Mno.06 spin glass sample 



is rapidly quenched to a temperature T < T g in the 
presence of a small magnetic field. The field is cut 
at time t = t w , and the magnetization decay is then 
recorded [III Si] for t > t w . 

The parameters shown in Table (1) are used empiri- 
cally determine, on the basis of Eq. [5J the SA exponent 
/i(T). The small deviation of the experimental data from 
the pure aging form given in EqJ3] implies, of course, a 
small systematic error. 

Of special importance are the values of the dominant 
exponent A2 which are listed in the last column of the ta- 
ble: According to Eq. 1101 the quantity A2 + x determines 
the time scale for the onset of EOA behavior. 

Figure [2] is scaling plot of the rate of (de-) magneti- 
zation multiplied by the system age, versus the scaling 
variable t/t w . At the aging temperature T = 0.83 the 
empirical sub-aging exponent /1 is very close to one 0, H| 
and the data conform reasonably well to a PA ansatz. 
The approximation is worse, but still usable, at other ag- 
ing temperatures and all data can be fit reasonably well 
by the full aging formula [3] 
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